Scaling of the conductance in gold nanotubes 
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A new form of gold nanobridges has been recently observed in ultrahigh-vacuum experiments, 
where the gold atoms rearrange to build helical nanotubes, akin in some respects to carbon nan- 
otubes. The good reproducibility of these wires and their unexpected stability will allow for conduc- 
tance measurements and make them promising candidates for future applications . We present here 
a study of the transport properties of these nanotubes in order to understand the role of chirality 
and of the different orbitals in quantum transport observables. The conductance per atomic row 
shows a light decreasing trend as the diameter grows, which is also shown through an analytical 
formula based on a one-orbital model. 
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Gold is known to be a good conductor for ages but its 
use in nanoelectronics delivers unexpected behaviors dif- 
fering from its bulk properties. Recent experiments^^ 
indicate that gold nanowircs exhibit fascinating ordered 
structures that resemble those of the by now well-known 
carbon nanotubes (CNTs)i^ 

High resolution microscope images show that these wires 
consist of coaxial nanotubes of helical atom rows coiled 
around the wire axis. These wires present therefore dif- 
ferent possible helicities as CNTs do, depending on the 
angle in which the atom rows rotate around the tube 
axis. Experimental evidence so far seems to support the 
fact that small chiral angles are preferred by the fab- 
ricated nanotubes. The multi-shell nanotubes keep the 
difference in the number of atom rows between adjacent 
shells constant, leading to a kind of magic shell-closing 
number which is always seven. This fact can be under- 
stood as a new inner shell arises when the radii differ 
in about a = 2.88 A, the neighboring distance of bulk 
gold, which implies a perimeter difference (2na ~ 7a) 
corresponding to seven atom rows nearly parallel to each 
other. So single- walled nanotubes should exist with six or 
less atom rows, and actually one single- walled tube with 
five atom rows has been observed^ With seven or more 
atom rows a double- walled tube should be seen, which be- 
comes a triple- walled nanotube when more than fourteen 
rows are present, and so on. The good reproducibility 
of these wires, and their stability^*SiL& which is better 
for those with outer tubes with an odd number of atom 
rows, makes them promising candidates for future appli- 
cations. Ballistic conductance measurements will lead to 
a confirmation of the proposed structural model through 
its derived electronic structure. Astonishingly enough 
several works hint at the fact that the conductance per 
atom row might decrease with the number of circumfer- 
ential atomsAifi In this Letter, we provide a theoretical 
understanding of this phenomenon by means of both nu- 
merical and analytical calculations. 

Determining the structure of the nanotubes is funda- 
mental in order to understand their physical properties, 



such as conductance quantization. Experimental evi- 
dence points out that a structural transition from the 
thicker, internally crystalline wires into these thinner reg- 
ular but noncrystalline ones takes place at a critical ra- 
dius below 2 nm. Lattice spacing has been measured to 
be almost 2.88 A, the neighboring distance of gold, for all 
these wires with diameters between 0.5 nm and 1.5 nm. 
We therefore take this distance as the interatomic spacing 
of the triangular lattice, as we can think of these cylin- 
drical tubes as rolled- up lattice planes of fee Au (111), as 
depicted in Fig. ^ Atoms in successive shells appear to 
strain a bit to maintain commensurability between inner 
and outer shells. Ab-initio calculations show though that 
this shear strain should have just a very small effect on 




FIG. 1: Some gold nanotubes are shown as well as the 2D 
triangular network, from which we can build these models. 
The coordinate system used is sketched on the lattice, being 
ai and a.2 the basis vectors, and C the chiral vector specifying 
the nanotube and corresponding to its circumference. The 
highlighted atoms on the lattice are those defining the region 
of inequivalent chiral vectors, which spans an angle of 30 ° 
as can be easily understood analyzing the symmetry of the 
lattice. 



Typeset by REVTeX 



2 



> 



a, o - 




DOS (arb. u.) 



G (2e /h) 



FIG. 2: Onedimensional band structure, density of states 
and conductance for the (4, 3) gold nanotube. For this tube 
the translational vector is T = |T| = 3.03 nm. The energies 
are shifted to set the origin at the Fermi energy, as in the next 
figures. 
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FIG. 3: Conductance of AuNTs with seven atom rows but 
different chiralities. Red lines indicate results including s or- 
bitals only, whereas the black ones are showing the calcula- 
tions with all nine outer orbitals. 



the electronic structure m^ 

Despite the strong intershell interactions, these struc- 
tures resemble greatly those of CNTs, where the hon- 
eycomb network of carbon atoms is replaced by the tri- 
angular lattice of gold atoms. These are actually comple- 
mentary networks, and the gold nanotube structure can 
be obtained from the one of CNTs by putting gold atoms 
at the center of honeycomb framework. 
We use the indices (n, m) to classify the gold nanotubes 
(AuNTs) , where n and m are integers defining the chiral 
vector, which is wrapped to form the tube. The chiral 
vector is then (nai + maa), perpendicular to the tube 
axis, where the vectors ai and a.2 span the 2D unit cell. 
The number of atom rows coiling around the tube axis 
is then n + m while the nanotube translational vector 
reads T(n,m) = ((2m + n)/dn) ai — ((2n + m)/dR) a%, 
being du the greatest common divisor of (2m + n) and 
(2n + in) , as T should be the smallest lattice vector in its 
direction. Due to the symmetry of the triangular lattice 
the chiral vectors of all inequivalcnt tubes are comprised 
in a 30 angle, that is, in a one- twelfth irreducible wedge 
of the Bravais lattice. We can then compromised and use 
only indices with n > and < m < n. The tubes of 
the form (n, 0) and (n, n) are achiral, presenting thus no 
handedness. Only wires with an even number of strands 
may have the structure (n, n)H 

For a fixed n, the radius is minimized for higher values 
of m (tension decreases with shrinking radius), so as to 
fulfill simultaneously tightest external packing and mini- 
mal wire radius £ The largest m for odd number of atom 
rows is n— 1, leading always to a finite chirality and thus 
a helical structure in these cases. The outer tubes with 
an even number of atom rows seem to take non-hclical 



structures with shorter periods. The helicity of the tubes 
is observed experimentally through a wavy modulation of 
the STM images^ 

We study the transport properties of several of these 
AuNTs using a multiorbital tight-binding Hamiltonian 
and Green function techniques. In particular we adopt 
a Slater-Koster-typei^ tight-binding approach, with the 
parameters taken from the nonorthogonal parametriza- 
tion of Papaconstantopoulos and co-workers mm. Our ap- 
proach remains of microscopic nature since the symme- 
tries of the atomic s,p, and d valence orbitals are taken 
into account. Ab initio studies including shell effects 
could be needed for greater quantitative accuracy. 
As a result we analyze the differences in the transport 
calculations that raise from the inclusion of all outer s, 
p and d orbitals instead of just the simplified model of 
s orbitals, which is a tempting approximation for these 
kind of systems and has been used in other works m In 
order to compare results, we will thus use both a one- 
orbital Hamiltonian describing just the s orbitals, and a 
nine-orbital Hamiltonian, containing s,p and d orbitals 
including all conduction electrons. In doing that we re- 
strict our calculations to nearest-neighbors interactions. 
The Hamiltonian describing our systems can be written 
as 



ft = X! H i<x,30 C i<x C j0 



(1) 



where the indices i,j run over the atom sites and the 
indices a, f3 mark the different orbitals. The transfer in- 
tegrals Hiajp are the on-site energies or hopping param- 
eters, dependent on whether i equals j or not. 
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FIG. 4: Conductance of AuNTs with different radii r. Red 
lines indicate results including s orbitals only, whereas the 
black ones are showing the calculations with all nine outer 
orbitals. 



The conductance is calculated from the transmission val- 
ues applying the Landauer formula, using Green func- 
tion techniques, 15 ^ In particular, we derive the elas- 
tic linear response conductance via the Fisher-Lee for- 
mula for the quantum mechanical transmission^: G = 

^Tr {r L g r R gi} , where r L/R = i (e l/r - S[ /R ) be- 

ing E L / R the self energy of the left or right lead respec- 
tively, and where Q is the Green function of the central 
region dressed by the electrodes. 

We consider semi-infinite single-walled nanotubes 
(SWNTs) as leads. Considering the leads as bulk gold 
will yield an increased imaginary part of the lead self 
energies, which will slightly decrease the conductance 
values. For multi-walled nanotubes (MWNTs) it has 
been suggested that the conductance G is mainly 
determined by the outer shell which is also responsible 
for the structural stability^ A slight reduction of the 
conductance is expected from the finite size effects. 
The calculated density of states and the conductance for 
the gold nanotube (4, 3) are shown in Fig. El This tube is 
one of the simplest possible candidates with seven atom 
rows around its circumference (Fig. The colored ar- 
eas show the contribution of the different orbitals to the 
total density of states, indicating a dominance of s or- 
bitals around the Fermi level. Considering only s orbitals 
around the Fermi energy is thus a good approximation, 
but when applying a bias it can be critical to include all 
s,p,d outer orbitals, as they play an important role in 
the opening of new conduction channels. 
A comparison of the conductance of nanotubes of the 
same diameter but different chiralitics, as shown in Fig. [31 
suggests that no significant difference is expected in their 
conductance values at the Fermi energy. For applied bias 




FIG. 5: Schematic illustration of the reciprocal lattice of 
the 2D hexagonal Bravais lattice of Au (111) where the first 
Brillouin zone is highlighted through a darker background. 
The k- vectors lying at the Fermi surface give rise to a nearly 
perfect circle in this 2D structure. The black cutting lines 
represent the Brillouin zone of a nanotube in an extended 
zone scheme which allows the use of the dispersion relation 
calculated for the 2D gold layer to get the bands of a folded 
nanotube. In this example we show the wave vectors giving 
the bands in a (6,0) AuNT, a good candidate for this visu- 
alization as it has a small number of bands. These lines are 
plotted as solid bands inside the first Brillouin zone, and col- 
ored when folded back into it through the reciprocal lattice 
vectors of the gold layer. In the inset the conductance per 
atom row at the Fermi energy based on the analytical model 
accounting only for s-orbitals described in this Letter is shown 
for all possible (n,m) NTs, reproducing exactly the values of 
the numerical calculations. 

i 

voltages the current will show the difference in energies 
at which new channels start taking part in transport. 
In Fig. 0] the conductance is plotted for nanotubes of 
different radii, chosen as to have the smallest chirality 
possible for a given radius. The number of atom rows 
around the nanotube axis in each of them is 6, 7, 9, 11, 13 
and 14 respectively. We see that, though the conductance 
increases for thicker nanotubes, this is not the general 
trend if divided by the number of atom rows. We can 
observe again how the expansion of the orbital basis to 
include all conduction electrons changes the energies at 
which new channels open, being critical to obtain correct 
current values. 

By analyzing the Brillouin zone of the two-dimensional 
gold lattice taking into account only s orbitals, one 
can easily demonstrate after studying its dispersion rela- 
tion that the Fermi surface is with great accuracy ap- 
proximated by a circle with a radius close to kp = 
2y / 7r '/ '(v^3 a) , where a = 2.88 A. This simple model allows 
us to have an analytical approach to the conductance, as 
the number of energy bands crossing the Fermi surface, 
since the Brillouin zone of the gold nanotubes consist of 
parallel line segments due to the quantization of the wave 
vector along the circumferential direction. For simplicity 



4 




infinite system 
with artificial barriers 
with wide band leads 



-1.5 



-0.75 





E F (eV) 



0.75 



1.5 



FIG. 6: Conductance of the (4,3)-AuNT for a perfect infinite 
NT (black line) , for a finite tube with barriers in the connect- 
ing regions to the electrodes (where the constant a introduced 
in the text equals 0.6) and for a finite tube like in the previ- 
ous case but in the wide band limit approximation taken at 
the Fermi energy. The colored lines show oscillations due to 
scattering at the ends of the finite size tube. The thicker lines 
correspond to room temperature conductance^- 



we can restrict this counting to the first Brillouin zone 
by folding back to this area the segments lying outside 
it, which will elongate the lines inside it, as pictured in 
Fig. The conductance can then be written as 



G = 



2e 2 



2Int 




(2) 



and is only dependent of the nanotube indices and not 
on hopping parameters or on-site energies. This mod- 
els yields a perfect matching with the numerical calcula- 
tions described before. The result of this approximation 
of the conductance at the Fermi level is plotted in the 
inset of Fig.[S]for nanotubes of different radii and chirali- 
ties, presenting the conductance per atom row. The solid 
lines in this figure put an accent on the results of AuNTs 
with the smallest helicity as characterized in Fig. We 
can observe how the envelope of all these points slowly 
decreases with increasing number of coiling atom rows, 
or likewise with increasing diameter. The conductance 
values will nevertheless be reduced by considering more 
realistic gold leads. This trend has been also found in a 
first-principles one-orbital calculation^ 



Gold nanotubes have been observed as short (L ~ 4 nm) 
conductors between bulk gold electrodes. This calls for 
an extension of our theory to account for finite size ef- 
fects. We did this by parametrically introducing into the 
nanotube two tunneling barriers at a distance L. Fig. El 
shows the conductance between two homologous semi- 
infinite electrodes which are taken into account through 
modified self energies. This modification is introduced 
by a multiplicative factor a, which simulates the bar- 
riers. The finite size effects can be equally made ap- 
parent by considering wide band leads, i.e. Hwb{E) = 
— iImS(i? = E-p). As we can see in Fig. El we obtain the 
same kind of oscillating behavior with both approxima- 
tions. This shows that the results of our model provide 
an upper limit for the conductance. 

In conclusion, we investigated the electronic properties of 
several gold nanotubes. We built the Hamiltonians of the 
systems by applying a tight-binding model for s orbitals 
as well as for all outer orbitals. All AuNTs are equally 
metallic independently of their chirality or diameter in 
contrast to the results obtained for their carbon counter- 
parts^ but expected from gold electronic nature. In the 
case of gold the quantization of electron waves along the 
circumference of the tube results in a Brillouin zone com- 
posed of line segments that always cut the Fermi surface 
of the Au [111] layer due to the continuity of this surface. 
The s-orbital calculations demonstrate that a one-orbital 
approach to the problem of gold nanotubes is not a bad 
approximation, but should be avoided when applying a 
finite bias voltage to the system. 

An analytical formula for the conductance of the s-orbital 
calculations, matching perfectly the numerical results, 
shows that there is a slight decrease of the conductance 
per atom row as the radius increases. 
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